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"cs" ! Abstract 

^^ ' Within the framework of perturbative QCD approach, we study the charmless two-body de- 

^\j ■ cays B — ?• ai(1260)iir*,6i(1235)-fC*. Using the decays constants and the hght-cone distribution 

Vh ' amphtudes for these mesons derived from the QCD sum rule method, we find the following 

results: (a) Our predictions for the branching ratios are consistent well with the QCDF results 

within errors, but much larger than the naive factorization approach calculation values. BarBar 

C""^ ■ has searched the decays B — t- a^K*^,biK* and set the upper limits for their branching ratios, 

which are (much) lower than our predictions. The current data seem to imply that penguin 

^ . annihilation is small in these penguin-dominated decays. Certainly, if the factorizable anni- 

Q^i hilation diagram contributions are turned off, the branching ratios will fall dramatically and 

O \ many of them reduced by more than half. It needs further accurate experiments to clarify the 

1^ ' differences between the theoretical predictions and the present upper limits, (b) We predict 

that the similar anomalous polarizations occured in decays B — )• (j)K* also happen in decays 

B — )• aiK* , while not happen in decays B — )• biK*. Here still the contributions from the anni- 

^ , hilation diagrams play an important role to explain the lager transverse polarizations in decays 

^^ [ B — )• aiK*, while are not sensitive to the polarizations in decays B — )■ biK*. (c) Our predic- 

Q"^ ' tions for the direct CP-asymmetries agree well with the QCDF results within errors. The decays 

'^ ■ B^ —7- b^ K*^ , B^ —7- b^K*^ have larger direct CP-asymmetries, which could be measured by 

cn . the present LHCb experiments. 
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I. INTRODUCTION 

In general, the mesons are classified in J^^ multiplets. There are two types of 
orbitally excited axial- vector mesons, namely l"^"*" and 1^ . The former includes 
ai(1260), /i(1285), /i(1420) and Kia, which compose the ^Pi-nonet, and the latter in- 
cludes fei(1235), /ii(1170), /ii(1380) and Kib, which compose the ^Pi-nonet. There is an 
important character for these axial- vector mesons except ai(1260) and 61 (1235), that is 
each different flavor state can mix with one another, which comes from the other nonet 
meson or the same nonet one. There is not mix between ai(1260) and 61 (1235) because 
of the opposite C-parities. They do not also mix with others. So compared with other 
axial-vector mesons, these two mesons should have less uncertainties about their inner 
structures. 

Like decay modes B -> VV , the charmless decays B -^ ai(1260)i^*, 61(1235)/^* 
also have three polarization states and so are expected to have rich physics. In many 
B — )■ VV decays, the informations on branching ratios and polarization fractions among 
various helicity amplitudes have been studied by many authors [11-14] • Through po- 
larization studies, some underling helicity structures of the decay mechanism are pro- 
claimed. They find that the polarization fractions follow the naive counting rule, that 
is /l ~ 1 — 0(my/m^),/|| ~ /± ~ 0{my/m\). But if the contributions from the fac- 
torizable emission amplitudes are suppressed for some decay modes, this counting rule 
might be modified in some extent even dramatically by other contributions. For exam- 
ple, many anomalous longitudinal polarization fractions in the decays B — t- pK*,(f)K* 
have been measured by experiments, which are about 50% |5[, except that of the decay 
B~ — )■ K*~p^ with large value (96^^g)% J5[ (the newer measurement is (90 ± 20)%) Q. 
The branching ratios of the decays B~ — )■ p~K*~ and B — )■ (pK* have been measured 
by experiments which are near 10.0 x 10~^ |5|, and that of the decay 5° — t- K*^p^ is 
smaller and about (3.4 ± 1.0) x 10~^ [5]. Whether the similar results also occurs in the 
B — )■ ai{1260) K* , bi{1235) K* decay modes is worth researching by theories and experi- 
ments. We know that ai(1260) has some similar behaves with the vector meson, so one 
can expect that there should exist some similar characters in the branching ratios and the 
polarization fractions between decays B — )■ ai(1260)i^* and B -^ pK*, where ai(1260) is 
replaced by its scalar partner p. While it is not the case for 61 (1235) because of its differ- 
ent characters in decay constant and light-cone distribution amplitude (LCDA) compared 
with those of ai(1260). For example, the longitude decay constant is very small for the 
charged 61 (1235) states and vanishes under the SU(3) limit. It is zero for the neutral 
6°(1235) state. While the transverse decay constant of ai(1260) vanishes under the SU(3) 
limit. In the isospin limit, the chiral-odd (-even) LCDAs of meson 6i(1235) are symmetric 
(antisymmetric) under the exchange of quark and anti-quark momentum fractions. It is 
just contrary to the symmetric behavior for ai(1260). In view of these differences, one 
can expect that there should exist very different results between B — )■ ai{1260)K* and 
B -^ bi{1235)K*. On the theoretical side, B -^ ai{1260)K*, bi{1235)K* decays have been 
studied by Cheng and Yang in Ref. |7| where the branching ratios are very different with 
those calculated by the naive factorization approach [8|. To clarify such large differences 
is another motivation of this work. On the experimental side, only the upper limits for 
some of the considered decays can be available |9|, [lO| . 



In the following, ai(1260) and 61 (1235) are denoted as ai and 61 in some places for 
convenience. The layout of this paper is as follows. In SecJTTl we analyze these decay 
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FIG. 1: Diagrams contributing to the decay B^ — >■ 05-^**^. 



channels using the PQCD approach. The numerical results and the discussions are given 
in Sec. Illli The conclusions are presented in the final part. 



II. THE PERTURBATIVE QCD CALCULATION 

The PQCD approach has been pro ved been an effective theory to handle hadronic 
B decays in many works [2|, y, [uJ, II4]. Because of taking into account the transverse 
momentum of the valence quarks in the hadrons, one will encounter double logarithm 
divergences when the soft and the collinear momenta overlap. Fortunately, these large 



double logarithm can be re-summed into the Sudakov factor [13|]. There are also another 



type of double logarithms which arise from the loop corrections to the weak decay vertex. 
These double logarithms can also be re-summed and resulted in the threshold factor. This 
factor decreases faster than any other power of the momentum fraction in the threshold 
region, which removes the endpoint singularity. This factor is often parameterized into 



a simple form which is independent on channels, twists and fiavors [ij]. Certainly, when 



the higher order diagrams only suffer from soft or collinear infrared divergence, it is ease 



to cure by using the eikonal approximation [15|. Controlling these kinds of divergences 



reasonably makes the PQCD approach more self-consistent. 

Here we take the decay B^ — )■ a^K*^ as an example, whose part of diagrams are shown 
in Figure 1. These eight Feynman diagrams belong to the condition of K*^ meson being 
at the emission position. For the factorizable and nonfactorizable emission diagrams [U^a) , 



[U^b) and [11(c), [11(d), if one exchanges the positions of a? and K*^, the corresponding dia- 
grams can also contribute to the decay. All of these single hard gluon exchange diagrams 
contain all of the leading order contributions to B^ — )■ a^K*^ in the PQCD approach. 
The amplitudes of diagrams [Tj^a) and [T]^b) are denoted as Fgai- If d{d) quark in meson 
qO ^^*o-) Qf each diagram in Fig.l is replaced with u{u) quark, we can get the Feynman 
diagrams for the decay B^ — )■ afK*~. It is ease to think that the decay B^ — )• afK*~ 
would receive larger tree operator contributions than the decay B^ — )■ aiK*^. If d quark 
in meson a^ (5°) of each diagram in Fig.l is replaced with u, we can get the Feyn- 
man diagrams for the decay B^ — )■ aiK*^. Compared with the decay B^ — )■ aiK*^, 
B^ — > a^K*^ and B~ — >■ a^K*~ have extra contributions from the the factorizable and 
nonfactorizable emission diagrams. It is similar to the decay modes B — )■ K*p |2|], both 
longitudinal and transverse polarizations can contribute to the decay width. So we can 
get three kinds of polarization amplitudes Ml (longitudinal) and Mn^t (transverse) by 
calculating these diagrams. Because of the aforementioned distribution amplitudes of the 
axial-vectors having the same formats as those of the vectors except a factor, so the for- 
mulas of here considered decays can be obtained from the ones of 5 — )■ VV decays by 
some replacements. Certainly, there also exists a difference: if the emitted meson is bi for 
the factorizable emission diagrams, the amplitude F^b-^ contributed by the {V — A)(y ±A) 
operators would be zero due to the vanishing decay constant fh^ . Here for our considered 
decays, the amplitude for the factorizable annihilation diagrams Faai (if cli meson is re- 
placed with 61, the amplitude is denoted as -FafeJ plays an important role. To calculate 
accurately them will have an important effect on the branching ratios and the polarization 
fractions. 



III. NUMERICAL RESULTS AND DISCUSSIONS 

For the wave function of the heavy B meson, we take 

^b{x, b) = -T^i^B + mB)750B(x, b). (1) 

Here only the contribution of Lorentz structure (l)B{x,b) is taken into account, since the 
contribution of the second Lorentz structure (pB is numerically small [l6| and has been 
neglected. For the distribution amplitude (pB^x^b) in Eq.([T|), we adopt the following 
model: 



M^x^ 1 



Mx,b) = NBx\l-x)'exp[~^^ - -{Ubbn (2) 



where lo^ is a free parameter, and taken as lo^ = 0.4 ± 0.04 Gev in numerical calculations, 
and Nb = 91.745 is the normalization factor for Ub = 0.4. 

In these decays, both the longitudinal and the transverse polarizations are involved for 
each final meson. Their decay constants, Gegenbauer moments and wave functions are 



the same with those in Ref. 17 



We use the following input parameters in the numerical calculations 18|, ll9| : 



/b = 190MeV, Mb = 5.28GeV, Mw = 80.41G'eV, 

tb± = 1.638 X 10~^^s,rBO = 1.525 x IQ-^^s, 

\Vus\ = 0.2252, \Vts\= 38.7 x 10-^ 7 = (67.2 ± 3.9)°, 

iKfel = 3.89 X 10-^ \Vtb\ = 0.88. 

In the B-rest frame, the decay rates of 5 — )■ ai{bi)K*, can be written as 
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where A^'^ is the total decay amplitude of each considered decay. The subscript a is 
the helicity states of the two final mesons with one longitudinal component and two 
transverse ones. The decay amplitude can be decomposed into three scalar amplitudes 
a, b, c according to 
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where M2 and M3 are the masses of the two final mesons ai(6i) and K*, respectively. The 
amplitudes A^^, A^at, A^t can be expressed as 



Ml = a tl{L) ■ tl{L) + 
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-el{L) . P,el{L) . P2. 
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(9) 



We can use the amplitudes with different Lorentz structures to define the helicity ampli- 
tudes, one longitudinal amplitudes Hq and two transverse amplitudes H±: 



Ho = MIMl, H± = MIMn T M2M3Vr^ - IMt, (10) 

where the ratio r = P2-P3/{M2M^). After the helicity summation, we can get the relation 

J2 M''^M'' = \Ml\^ + 2{\Mn\^ + \Mt\^) = \Ho\^ + \H+\^ + \H_\^. (11) 



=L,N,T 



Certainly another equivalent set of helicity amplitudes are often used, that is 

Ao = -MIMl, 

All = V2MIMn, 

A^_ = M2M3v/2(r2 - 1)Mt. 

Using this set of helicity amplitudes, we can define three polarization fractions /o, 

\ A u |2 



/o,| 



Uo|2+U|||2+Uj 



(12) 



(13) 



The matrix element A4j of the operators in the weak Hamihtonian can be calculated 
by using PQCD approach, which are written as 

= V^,V:j,{l + z,e^^^^'^^), (14) 

where j = L,N,T and 7 is the Cabibbo-Kobayashi-Maskawa weak phase angle, defined 
via 7 = arg[— y''y^ ]. 5j is the relative strong phase between the tree and the penguin 
amplitudes, which are denoted as "T," and "-Pj", respectively. The term Zj describes the 
ratio of penguin to tree contributions and is defined as 

,.- Yl^ ^ (15) 

In the same way, it is easy to write decay amplitude Aij for the corresponding conjugated 
decay mode: 

M, = V:,VusT, - V^lVtsP, 

= K;KsT,(l + z,e^(-^+^^)). (16) 

So the CP-averaged branching ratio for each considered decay is defined as 

B=i\M,\' + \M,\')/2 = \VubV:f [Tlil + 2zLCos^cos5L + zl) 

+2 ^ Tf{l + 2zj cos 7 cos Sj + z^) 



=N,T 



:i7) 



Like the decays of B to two vector mesons, there are also 3 types of helicity amplitudes, so 
corresponding to 3 types of Zj and 6j, respectively. It is easy to see that the dependence 
of decay width on 6 and 7 is more complicated compared with that for the decays of B 
to pseudoscalar mesons. 

Using the input parameters and the wave functions as specified in this section, it 
is easy to get the branching ratios for the considered decays which are listed in Table 
m where the first error comes from the uncertainty in the B meson shape parameter 
Ub = 0.40 ± 0.04 GeV, the second error is induced by the hard scale-dependent varying 
from ^QCD ~ 0-25 ± 0.05, and the last one is from the threshold resummation parameter 
c varying from 0.3 to 0.4. In FigJ21 we also show the Cabibbo-Kobayashi-Maskawa angle 
7 dependence of the branching ratios of decays B — > aiK* and B — )■ biK*. It is easy to 
see that the branching ratios of decays B^ — )■ ai{b^)K*~ are more sensitive to the angle 
7 compared with those of other decays. Several remarks on the numerical results are in 
order: 

• Just like the decays 5° -^ afp^^aip^), the decays 5° — > a^ K*^ {a^K*'^) have very 
different tree contributions (shown in Table |Tl]), for one is color allowed and the 
other is color suppressed. While the former proceed via 6 — )■ c? transition, where 
the tree operator contributions are CKM allowed, so the branching ratio of a^p~ 
is much larger than that of the decay a^p^ |7|. The latter proceed via 6 — )■ s 
transition, where the tree operator contributions are CKM suppressed, so the decay 
modes afK*~ and a^K*'^ have the branching ratios in the same order, because their 
penguin contributions are near each other. 

6 



TABLE I: Branching ratios (in units of 10 ®) for the decays B — ^ ai(1260)-ff* and B — ^ 
6i(1235)i^*. In our results, the errors for these entries correspond to the uncertainties from lob, 

(5) 

the QCD scale A^^^ and the threshold resumniation parameter c, respectively. For comparison, 
we also listed the results predicted by QCDF approach [7] and the naive factorization approach 
The upper limits of the decays B — )• a7 K*^, biK* measured by BaBar are also listed. 
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FIG. 2: The dependence of the branching ratios on the Cabibbo-Kobayashi-Maskawa angle 7. 
The left (right) panel is for the decays B — )• ai{bi)K*. The dotted line represents the decays 
B- -^ a\{bl)K* 
the decays B^ - 



, the solid line represents the decays B^ — )• a1{b\)K*^ , 
a^{b^)K*^, the dot-dashed line is for the decays B^ — > 



the dashed line is for 
at{bt)K*-. 



In our prediction, the branching ratio of the decay B^ — )■ aiK*^ is larger than that 
of the decay B~ — )■ a\K*~ , it is mainly induced by the amplitudes of factorizable 
emission diagrams [It^a) and [U^b) Fea^ and F^k* (for the corresponding exchange 
diagrams) have contrary interference effects between these two decays : constructive 
for the decay a^K*^ , destructive for the decay a\K*~ . So the the decay B^ — )■ 
^0^*0 ];gcgives a larger real part for the penguin amplitudes. Though the decay 
B^ — )■ a\K*~ has much larger contributions from tree operators, which are CKM 



TABLE II: Polarization amplitudes of different diagrams for the decays B^ 
(xlO-^GeV^). 
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suppressed and can not change the branching ratio too much. 

In order to characterize the contribution from tree operators and symmetry breaking 
effects between B~ and B^ mesons, it is useful to define the two ratios: 



B{B- -^ a:[-K*^) r^o „ B{B- -^ h\K* 



r^o 



BiB'^^alK*-) tb- B{B^^h+K*-) r^- " ' 

If one neglects tree operators and electro-weak penguins, the ratios obey the limit 



R^ = l,i?2 = 0.5. 



(19) 



Here the two ratios -Ri and R2 are predicted as 1.02 and 0.55, respectively. It is 
to see that the two ratios predicted by QCDF are 0.98 and 0.52, respectively. If 
the future data for Ri have large deviation from our value, the contribution from 
electro- weak penguin operators might give an important affect, for the contribution 
from tree operators can not change the branching ratio of B^ — )■ afK*~ too much. 
If the future data for R2 have large deviation from our calculation value, some 
mechanism beyond factorization even from new physics might give an important 
affect, because the factorizaton formulae between B^ — )■ bfK*~ and B~ — )■ h\K*~ 
are exactly the same with the neutral h\ meson decay constant vanishing. 

Compared with other results: From Table [H One can find that our predictions 
are consistent well with the QCDF results within (large) theoretical errors, while 
in stark disagreement with the naive factorization approach. Our predictions are 
much larger than the results calculated by Calderon, Muiioz and Vera [8|, where 
the nonfactorizable effects are described by the effective number of colors N^^^ . For 
some decays, where the contributions from the emission diagrams are dominated or 
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TABLE III: Branching ratios (in units of 10~^) for the decays B —?- ai(1260)-ftr* and B —?- 
6i(1235)i^*. The label -Nonfac. Anni. and -Fac. Anni. mean the results with neglecting the 
nonfactorization and factorization annihilation diagrams, respectively. The label -Anni. means 
the results with neglecting all the annihilation type contributions. 





Full contributions 


-Nonfac. Anni. 


-Fac. Anni. 


-Anni. 


B^ -^ a+K*- 


9.9 


10.3 


5.2 


5.4 


BO ^ a?K*o 


7.1 


7.4 


5.0 


5.2 


B- -^ a{K*° 


10.8 


11.0 


6.7 


8.4 


B- -^ a\K*- 


4.8 


4.9 


2.1 


2.1 


50 -^ b+K*- 


18.0 


19.0 


8.5 


9.0 


50 ^ 50^*0 


9.6 


10.2 


4.5 


5.1 


B~ -^ b^K*^ 


23.0 


23.6 


10.9 


12.4 


B- -^ b\K*- 


10.6 


10.9 


4.8 


5.0 



the branching ratios have a strong dependence on the correlative form factors, the 
naive factorization approach can give a reasonable prediction, while for the decays, 
where the annihilation diagrams play an important role, this approach would expose 
some disadvantages. BarBar has been searched the decays B — )■ a^K*^ , biK* and set 
the upper limits on their branching ratios ranging from 3.3 to 8.0 x 10^^ at the 90% 
confidence level J9|,[l0|. Our predictions are (much) larger than the upper limits. It is 
noticed that these upper limits are obtained by assuming that B{ai — )■ 7r"'"7r~7r^) = 
B{af — )■ 7r°7r°7r='=) and B{af{bf) — )■ p°(a;)7r^) = 1. The current data seem to imply 
that penguin annihilation is small in these penguin-dominated decays. Certainly, 
if one neglects the contributions from the nonfactorizable annihilation diagrams, 
the branching ratios only have a tiny increase, while if the factorizable annihilation 
diagram contributions are neglected, the branching ratios will fall dramatically and 
many of them reduced by more than half (shown in Table IIIip . It needs further 
accurate experiments to clarify the differences between the theoretical predictions 
and the present upper limits. 

To resolve the B — )■ (pK* polarization puzzle surges a considerable amount of theoretical 
attentions: The author in Ref . |20| considered two class new-physics four-quark operators, 
namely (1 ± 75) ® (1 ± 75) , a(l ± 75) cr(l ± 75), which could account for the data. 
While the author in Ref. |2l| proposed to take a smaller B — )■ K* form factor and 
strengthen the penguin annihilation and nonfactorizable contributions, which could be 
achieved simultaneously by adopting the asymptotic models for the K* meson distribution 
amplitudes. Here we take the second suggestion and find that the polarization charactors 
for decays B — )■ aiK* and B — )• biK* are very different: the transverse polarization 
amplitudes have almost equal values with (even a little stronger than) the longitudinal 
polarization ones for the former, while the longitudinal polarization states are dominated 
for the latter. It seems that the similar anomalous polarizations occured in decays B — t- 
(pK*, pK* also happen in decays B — )• aiK* , while not happen in decays B — )• biK*. Here 
we find that the contributions from the annihilation diagrams are very important to the 
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TABLE IV: Longitudinal polarization fraction (/l) and two transverse polarization fractions 
(/||i f±) for the decays B — )• ai{12Q0)K* and B — )• bi{1235)K* . In our results, the uncertainties 

(5) 

come from ujb, the QCD scale Ag^^ and the threshold resummation parameter c. The results 
of /l predicted by the QCDF approach are also displayed in parentheses for comparison. 



/l(%) 



/||(%) 



f±{%) 



B^- 
B^ 
B- 
B- 



a{K 



48 
59 
a^/?*° 50 



1 



aO^*o 



a'^K*- 49 



Q+5. 1+7.4+4.9 
■^-4.7-8.0-4.9 
f^+4.7+7.7+4.3 
"-4.9-7.8-4.3 
q+5. 1+8.6+5.0 
•-'-4.9-9.9-5.0 
n+3. 3+6.2+4.7 
"-4.3-6.2-4.7 



-39 \ 
-29 J 



(37 

i23l:tl) 

48 ^ 
37 J 

(521:11) 



(37^ 



+2.5+3.8+2.6 
-2.8-4.1-2.6 
)+2.6+3. 8+2.2 
-2.5-3.8-2.2 
24 -|^+2.6+5. 0+2.5 



26.1] 

20.2^ 



25.5: 



-2.7-3.7-2.5 
1-2.3+0.0+2.4 
-2.5-2.5-2.4 



25.0 
20.2 
25.6 
25.5 



+2.2+3.8+2.3 
2.3-3.5-2.3 

+2.3+4.0+2.1 
2.2-3.5-2.1 

+2.3+5.0+2.5 
2.4-4.9-2.5 
2.0+3.2+2.2 
2.2-3.7-2.2 



B^- 

B" 
B 



bJK 



95 

h\K*^ 95 
6^^*o 96 
h\K*- 96 



Q+O.l+l. 1+0.0 
^-0.1-1.3-0.0 
4+0.1+1.0+0.1 
■^-0.1-1.4-0.1 
9+0.0+0.9+0.1 
^-0.0-1.7-0.1 
r-+0. 0+0.8+0.1 
'-'-0.1-1.3-0.1 



(82; 
(79: 
(79: 

(82: 



-18^ 
-41 J 
-2U 
-74 J 
-2U 
-74 J 
-16^ 
-26 J 



1.1 

0.9 
1.0 
0.7 



+0.2+0.4+0.2 
-0.0-0.2-0.2 

0.0+0.2+0.4 
-0.0-0.2-0.4 

0.0+0.3+0.3 
-0.0-0.3-0.3 

0.1+0.2+0.2 

0.0-0.1-0.2 



3.0 
3 

2 

2, 



+0.0+0.9+0.2 
0.1-0.7-0.2 
7+0.1+1.2+0.3 
'-0.1-0.8-0.3 
Q+0. 0+0.9+0. 2 
"-0. 0-0.6-0. 2 
Q+0.0+0. 9+0.3 
°-0.0-0.6-0.3 
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FIG. 3: Direct CP-violating asymmetry as a function of Cabibbo-Kobayashi-Maskawa angle 7. 
The dashed line is for the decay B~ — )■ h^K*^ , the solid line represents the decay B^ — t- b^K*^ , 
the dotted line represents the decay B^ — )■ h\K*~ ^ the dot-dashed line is for the decay B^ — )• 



bJK*' 



final polarization fractions for the decays B — )■ aiK*: If these contributions are neglected, 
the longitudinal polarization fraction of the decay B~ — )■ a\K*~ becomes 98.8%, those 
of i?° -)■ af K*~ , a'lK*'^ increase to about 90%, that of the decay B^ -)■ a^K*° changes 
from 50.3% to 70.0%. While the longitudinal polarizations of decays B^ — )• h^K*' and 
B^ — )■ bfK*~,biK*^ only have a very small decrease by neglecting the annihilation type 
contributions. That of the decay B~ — )■ b~[K*'^ has a larger amplitude reduction, changing 
from 96.2% to 86%. In a word, the longitudinal polarizations of decays B — )■ biK* are 
not sensitive to the annihilation type contributions. 
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Now we turn to the evaluations of the CP-violating asymmetries in PQCD approach. 
Here we only research the decays B — )■ biK*, where the transverse polarization fractions 
are very small and range from 3.8 to 5.2%. Using Eq.f ll4p and Eq. (TT6!) . one can get the 
expression for the direct CP-violating asymmetry: 

.dir ^ \M\^-\M\^ ^ 2zl sin a sin 6l .^q) 

^^ \M\^+\M\^ il + 2zLCosacos6L + zl)' ^ ' 

Here for our considered four decays, the contributions from the transverse polarizations 
are very small, so we neglected them in our calculations. Using the input parameters and 
the wave functions as specified in this section, one can find the PQCD predictions (in 
units of 10"^) for the direct CP-violating asymmetries of the considered decays: 

Atp^B"" -^ b+K*-) = 38.5tl-^^tflttl (21) 

A^UB- ^ b'.K*-) = 5A.3tr7tl',ttl, (22) 

A'MB' ^ 6?i^*°) = -18.7tlT-'or.ll (23) 

A^UB- ^ b^K*') = -OASt'oin^lt'oil (24) 

where the errors are induced by the uncertainties of B meson shape parameter Ub = 
0.4 ±0.04 , the hard scale-dependent varying from A^^^ = 0.25 ±0.05, and the threshold 
resummation parameter c varying from 0.3 to 0.4. It is particularly noteworthy that our 
predictions about the direct CP asymmetries of these decays are consistent well with the 



QCDF results |22 



A'SUB' ^ biK*-) = (44+^3)%, (25) 

A'SUB- ^ b',K*-) = (60+?3)%, (26) 

A'^UB'' -^ b'iK*'') = (-17+21)%, (27) 

A''jp{B- -^ b^K*"") = (2+°)%, (28) 

where the error comes from the parameters pA,H and arbitrary phases (I)a,h- These are 
phenomenological parameters to cure the endpoint divergences in the amplitudes for the 
annihilation and hard spectator scattering diagrams. Just because of cannot calculat- 
ing the annihilation and hard spectator scattering diagrams accurately, there are large 
theoretical errors for the QCD factorization predictions. 

IV. CONCLUSION 

In this paper, by using the decay constants and the light-cone distribution amplitudes 
derived from QCD sum-rule method, we research B — )• aiK*,biK* decays in PQCD 
factorization approach and find that 

• Our predictions for the branching ratios are consistent well with the QCDF results 
within errors, but much larger than the naive factorization approach calculation 
values. BarBar has searched the decays B — )■ a^K*'^, biK* and set the upper limits 
for their branching ratios, which are (much) lower than our predictions. The current 
data seem to imply that penguin annihilation is small in these penguin-dominated 
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decays. Certainly, if the factorizable annihilation diagram contributions are turned 
off, the branching ratios will fall dramatically and many of them reduced by more 
than half. It needs further accurate experiments to clarify the differences between 
the theoretical predictions and the present upper limits. 

we predict that the similar anomalous polarizations occured in decays B — )■ (pK* 
also happen in decays B — t- aiK*, while not happen in decays B — )■ hiK* . Here still 
the contributions from the annihilation diagrams play an important role to explain 
the lager transverse polarizations in decays B — ?■ aiK*, while are not sensitive to 
the polarizations in decays B — )■ hiK*. 

Our predictions for the direct CP-asymmetries agree well with the QCDF results 
within errors. The decays B^ — )■ bfK*~,B~ — )• h\K*~ have larger direct CP- 
asymmetries, which could be measured by the present LHCb experiments. 
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